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We find a solvable limit to the problem of the ID electron gas interacting with a lattice of Kondo 
scattering centers. In this limit, we present exact results for the problems of incommensurate filling, 
commensurate filling, impurity vacancy states, and the commensurate-incommensurate transition. 



PACS numbers: 

In this paper we describe an exactly solvable limit of 
a one dimensional electron gas interacting with a lattice 
of dynamic spin scattering centers. We are interested in 
this model both in its own right and for its possible ap- 
plicability to the understanding of Kondo insulators and 
heavy fermion materials. Among our most interesting re- 
sults, we find: 1) In addition to the single impurity Kondo 
scale, r, there is a coherence scale, A, which characterizes 
the temperature (or frequency) below which coherence 
between different impurities is established; this scale is 
also roughly equal to the gap, A s , in the spin excitation 
spectrum which appears at low temperatures regardless 
of whether the impurity lattice is commensurate, incom- 
mensurate, or even weakly disordered. Of course, A is 
an increasing function of impurity concentration, c, and 
vanishes in the c — ► limit. 2) There is a well defined 
crossover as a function of c from a dilute limit, in which 
A(c) <C r to a dense limit in which A(c) ^> T. In the 
dilute limit, single-impurity Kondo physics is apparent 
for temperatures in the range T > T > A, and coherence 
between impurities sets in only for A ~ T, while in the 
dense limit, coherence and Kondo quenching of the im- 
purity spins occur simultaneously, in much the same way 
as pairing and condensation occur simultaneously in BCS 
superconductors. 4) At the crossover concentration, c* , 
we find that A(c*) ~ T and c* ~ this confirms the 

existence of a physical length scale £k = vf/F associated 
with the single impurity Kondo problem. 5) We find that 
the coherent low-temperature state is characterized by 
the appearance of long-range order of a non-local order 
parameter, analagous to the hidden (or topological) order 
found in integer spin-chains, or the Girvin-MacDonald 
order parameter in the fracional quantum Hall effect. |^] 

In a discrete model, the free electron gas is described 
by a hopping Hamiltonian between lattice sites with lat- 
tice constant a. We consider the problem in which the 
electrons interact, via Kondo scattering, with a periodic 
array of localized dynamic spins with lattice constant b. 
The relative impurity concentration is given by c = a/b. 
We focus on the long-distance behavior of the electrons' 
correlation functions by taking the continuum limit of the 
electron lattice. At first we leave the interaction with the 



impurity lattice in its discrete form. The whole system 
is described by the Hamiltonian 



H = Hn + Hu 



(1) 



Here Hq is the kinetic energy of noninteracting electrons 
with spin, which, in the continuum limit is (setting h = 1) 



Hq = —wf 



dx 



*LA*H,«r - *LA*i.* . ( 2 ) 



where ^\ (x) with A = R, L create respectively a right 
and left moving electron with z component of spin a = 
±1/2 at position x. H\\ +H± is the coupling between the 
electron gas and the discrete array of anisotropic Kondo 
impurities (which we refer to as "impurity-spins" ) : 



J, A 



H 



j -l E [^/*A|(^-)*An(^)+H.c. 



J.A.A' 



(3) 



Here rj* represent the dynamics of the impurity-spin at 
position Rj. (ig satisfy canonical comutation relations 
I 7 "/' r i'] = *^Jii' e -' c Tj ). We will always work in the limit 
of large anisotropy, Jii 3> Jj_. 

Bosonization of the kondo lattice: We bosonize 
the Hamiltonian using the standard one dimensional re- 
lation between Bose and Fermi fields |TJ : 



zexp{i$x,a(x)} 



(4) 



where $a.<t = \^[&a(x) ± <frcr(x)] ± k*FX with "+" and 
"-" corresponding to A = Left and Right respectively, 
a (x) = f_ dx'H (T (x'), <j) a {x) and II ff (x) are canoni- 
cally conjugate Bose fields, so that [(j> a (x),Il <T >(z')] — 



i5(x 



To take advantage of the separation between 



charge and spin which occurs in one dimension, we ex- 
press the Hamiltonian in terms of a spin field, <p s (x) — 
[0f — </>i\/y/2, and a charge field, 4> c {x) — [<fo + 0j/v2 
and correspondingly defined momenta, Tl s and n c . We 
introduce distinct coupling constants which we label with 
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the superscripts "f" and "b" for the forward and back- 
scattering terms in the Hamiltonian. This distinction can 
be associated with the form factor of a spatially extended 
impurity. The result is easily seen to be: 

H = V -f J dx {[Til + (dAc) 2 ] + [H 2 a + (dAs) 2 }} (5) 

3 

+ — E T i sH^Mj)] sin[V2^4> c (j) + ZkpRj], 
j 

H ± =E r 7 fe " i ^°' ) {— cos[V2^^(i)] 
z — ' J ira 

3 

Jb 

+ cos[V2^> c (j) + 2k F Rj}} + H.c. 
ira 

We eliminate the explicit dependence on the 9 S field by 
making a unitary transformation [0,0] 



U = exp [-iVtor^TjO^Rj)]. 



(6) 



The resulting transformed Hamiltonian is H = IPHU 
H = H° + H c + Ajf a J?-J2 T R d Mx3), 



2 Jit 

#n = — E f ^~ l ) j 

air * — ' J 



x sin 



2nM)] sin[V27r0 c (j) + 2k F Rj 



Hj = 



J 



^E cos[V2^ a (j)] 
7ra * — ' 



Jb 

V ff cos[V2n(t) c (j) + 2k F Rj], 
na £ — ' 



(7) 



where ff = UrfU + , and AJ|| — — irvp. We see that 
the result of the unitary transformation is to adsorb a 
spin degree of freedom (the bosonic phase field 6 S ) into 
the impurity operator r + leaving behind a phase shift tt 
of the number field v27r</>s between impurity sites (giving 
rise to the staggered coefficient (— We interpret this 
as a description of the Kondo resonance in the bosoniza- 
tion language. Therefore we can interpret the unitary 
transformation as a mapping from a free-field to a strong 
coupling basis. In the new basis the Hamiltonian is so 
simple that much of the physics can be read directly from 
it. 

Incommensurate phase: The solvable limit [^]]|] 
Jjf = Vpir (i.e. AJii = 0) is analogous to the Toulouse 



limit of the single impurity Kondo problem. For incom- 
mensurate filling, 2k F Rj ^ n7r, we make the approxi- 
mation of neglecting all terms in the Hamiltonian which 
depend on 2kpRj. (We will return to this point, below, 
when we consider the commensurate-incommensurate 
transition). In this limit, 



Hinc — Hg 



irn — ' 



fx, 
3 



(-l)^cos[V2^(j)] +J ff c . (8) 



{ff} commute with the Hamiltonian, so we can regard 
them as c-numbers (ff — ±5), labeling eigenstates of 
the transformed Hamiltonian with static configuration 
of ff. Also, the charge and spin parts of the Hamilto- 
nian decouple. The charge excitations, <fi c , are just the 
gapless free field states of Hq. The spin Hamiltonian has 
a discrete sine-Gordon form. Due to the (— 1) J factors, 
the ground-state is one in which the transformed impu- 
rity spins ff are antiferromagnetically ordered and the 
4> s spectrum is gapped, E B (k) = ±y/ (vpk) 2 + A 2 .. With 
respect to the ground state of this reduced Hamiltonian, 
small perturbations in AJy and Jn are irrelevant. There- 
fore, results obtained in this limit are generic for at least 
a finite region of parameter space. 

Physically there are two important energy scales: The 
single impurity Kondo resonance energy scale T, (which 
in the Toulouse limit is T = Jj_/irupa), and the inter- 
impurity coherence scale A which is set by the spin gap, 
A s , calculated below. The interplay between the two 
energy scales is sensitive to the impurity concentration. 

The Dense Limit - Associated with the spin gap A s 
in the spectrum there is a corresponding fermionic cor- 
relation length £ s = Vp/A s . When £ s is much bigger 
than the characteristic distance between impurity sites, 
d = fl/c«( s , the discrete character of the imurity array 
can be ignored. Therefore, in this limit (high density), 
we can also take the continuum limit with respect to the 
impurity lattice. The spin Hamintonian now takes the 
regular sine- Gordon form 



H„ — Hn 



2na 2 



dx cos[/?</> s (x)], 



with P 



2tt. The resulting gap is M 



A fi 



Up 

a 



cJ]_/v F 



2/3 



(9) 



(10) 



Similarly, there is an energy gap A ~ A s for creation of 
kinks in the transformed impurity spin order. 

Thus, the system is dense if c exceeds the critical den- 
sity ci, at which £ s (ci) = a/ci; consequently, in this limit 
A s > r. Since A s also determines the temperature scale 
below which coherence sets in, this is analogous to the 
BCS limit of superconductivity, in that local coherence 
between the impurity spins and the conduction electrons, 
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and longer range coherence between the impurity spins 
occur at the same temperature. 

Dilute limit - When the distance between impurities 
is large, a/c a, the discrete nature of the impurity 
array cannot be ignored. For small enough impurity con- 
centrations c < C2, we can compute || the leading order 
dependence of A in powers of the distance between impu- 
rities; A = (ttvfc/Ao) [1 + 0(cvf/J±)]- The consistency 
of the dilute limit is determined by the simultaneous con- 
ditions, cvf/J± "C 1 and A -C T, of which the latter is 
the more restrictive by a factor Vf/J±. From this condi- 
tion we get C2 ~ (J±/itvf) 2 ■ Note that there is actually 
only one characteristic concentration, c\ ~ C2 ~ c*. In 
the low temperature coherent regime, T <C A, the be- 
havior of both dilute and dense systems are qualitatively 
indistinguishable; for simplicity, we will henceforth deal 
with systems in the dense limit. 

Order parameters - The most interesting physical 
quantities the correlation functions of the various pos- 
sible order parameter fields. (Since this is a one di- 
mensional system, we will never see a broken continu- 
ous symmetry, but we can interpret the emergence of 
fairly long range coherence in the fluctuations of a dom- 
inant collective field as the one dimensional analog of 
the phase transition to a broken symmetry state that 
we would expect in higher dimensions.) In the case 
of noninteracting electrons, density-density correlation 
functions decay as l/x 2 . Therefore, any of the various 
possible order parameters whose correlation functions 
Ci{x, x') =< Oi(x)Oi(x') > decay like x~ ai with on < 2, 
can be considered as having enhanced long range coher- 
ence. The enhanced order parameters are: 2kp charge 
density wave Ocdw — [^l t^R,t + *l i^ii.i]) compos- 
ite spin density wave O c sdw — [^jr ^ R^rf,-.] (where 
j = j(x) is the j-th site nearest to position x), singlet 
pairing S p 



O cT P = , * ',, i t 



[w L ,], and composite triplet pairing 

L lTjr x j]- Each exponent cti has contri- 
butions from the spin and charge correlation function 
exponents K* and K c as follows: o>cdw = K c + K*, 
asp = + K*, a cS DW = K c , and a cT p = 1/K C . 

Since the spin field is massive, the spin-fluctuations are 
frozen out, which results in K* = in these correla- 
tion functions, and causing exponential decay of others. 
Moreover, deep in the incommensurate phase the charge 
fields are effectively free, so K c = 1. Note that the uni- 
tary transformation has eliminated the spin field depen- 
dence of the composite order parameters. 

All the above order parameters have staggered order 
Ci(x,x') ~ (— >\x — x'\~ ai . The impurity corre- 

lation function < iftfi > decays exponentialy. How- 
ever, the transformed impurity spin correlation function, 
< TjTj, >~ const. , exhibits long-range order at T = 0; 
this is the non-local order parameter which characterizes 
the coherent state. || 

The composite triplet pairing is especially interesting. 



Order parameters of this type were previously found in 
the solution of the single impurity 2-channel Kondo prob- 
lem 0], but not in the single channel Kondo problem. In 
going to the incommensurate limit we have neglected all 
the backscattering terms. If we consider the Kondo in- 
teraction without backscattering then it can be written 
as 



H impur = J f J2 T R ■ \M a (RjW^(Rj) 



(11) 



where i=R,L. The impurity scatters the right and left go- 
ing electrons as if they were two independent channels. 
Of course, for the single impurity problem, J b would lift 
this artificial channel degeneracy, and hence is a relevant 
pertubation. Even if J* <C J h , the single impurity is ulti- 
mately governed by the single-channel fixed point. How- 
ever, for the array, the existence of the spin gap stabilizes 
the two-channel behavior of the model with J b = 0, ren- 
dering small J b an irrelevant perturbation. We would like 
to note that a composite triplet pairing order parameter 
of a different kind was found in a mean field theory of 
isotropic the 3D Kondo lattice The relation between 
the two results is not clear at the present. 

Commensurate filling: In the case where fop is 
commensurate with the impurity lattice, the back scat- 
tering terms cannot be ignored. In particular we will now 
discuss the case 2kpR^ = jit (usually called half- filling). 
The transformed Hamiltonian is then 

H = H s + — V ff(-l)> J[ cos[^ s (j)] 
ira * — ' 

3 

+ H o + — E f /(-i) J J i cos[VW) C (j)] 



2Jf 



if sm[V2^Mj)} sin[VW c (j)] (12) 



0, ff are good quantum numbers. 



In the limit j| - ... , ( 

The spin and charge decouple into two independent sine- 
Gordon Hamiltonians, which are minimized when the <p 
fields take the values 2im. Consequently, both spin and 
charge excitation spectra are gapped, and ff are ordered 
antiferromagnetically. The quantum numbers of the ex- 
citations are easily determined in the bosonization for- 
malism. Solitons in the 4> s field are ^/2tt phase slips, 
amounting to a state with spin S = 1 and charge Q = 0. 
By the same analysis, a <fi c soliton carries 5 = and 
Q = 2e. A kink (or domain wall) in the ff order induces 
compensating ^V2n phase slips in cf> c and (f> s , resulting 
in a bound state with charge Q=l and spin S=l/2. As 
could be expected for commensurate filling, the only or- 
der parameter which does not decay exponentially is the 
2k F CDW order; Ocdw = [^*KT + *L*fl±]- At long 
distances < Ocdw(x)Ocdw( x ') >^ constant, since all 
fields <p s ,(f> e ,fj are gapped. Deviations Jy ^ can be 
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treated perturbatively giving rise to short range hopping 
dynamics of fj kinks. 

Note that there is a second possible stable phase when 
Jy 7^ in which r z orders and the </> fields take the val- 
ues 2-nn +1/2. This second stable phase is analagous to 
the dimerized commensurate phase of a half-filled Peierls 
insulator, || such as polyacetylenc; its implications will 
be explored elsewhere. 

Vacancy states: What happens when we introduce 
imperfections by deleting a single impurity from the pe- 
riodic array? The perfect array described by the Hamil- 
tonian in equation (12) had a staggered factor (— l)- 7 in 
front of both cos[\/27r</> s (j)] and cos[-\/27r</> c (j)] terms, but 
the origin of this factor is different for each term. The 
(— iy factor in the </> s part of the Hamiltonian was a 
result of the unitary transformation; it "counts" impuri- 
ties, but is insensitive to the distance between impurities. 
The (— 1) J in the <p c part came from the 2kpRj phase in 
coa[^/2n <j) c (j) + 2k F Rj], and therefore explicitly depends 
upon the distance between Kondo spins. Across a Kondo 
vacancy there will be a conflict between the impurity fj}. 
order induced by the spin fields and the charge field. This 
can be amended by either twisting the 4> s or the <p c fields 
across the vacancy site. Therefore, associated with a va- 
cancy there could be two alternative, mutually exclusive, 
bound states: Either a local </> c soliton of Q=l and S=0, 
or a local <fi s soliton with Q=0 and S=l/2. This is anal- 
gous to the solitons in polyacetylene. [|| 

Commensurate-Incommensurate Transition: 
For incommensurate wave number k F we can define in- 
commensurability q such that (k F +q)Rj = jir. Making a 
simple change of variables <p c — > <f> c — 2qxj \/2~k, the trans- 
formed incommensurate charge part of the Hamiltonian 
can be rewritten in the form 



so there is a significant range in which the soliton lattice 
exists. 

The gapless Goldstone modes, associated with the bro- 
ken symmetry soliton lattice ground state, are not triv- 
ially related to the free <j) c fields. Therefore, we can ex- 
pect correlation functions of operators like ^^ x ' (i.e. 
the order parameters) to have anomalous exponents. In- 
deed, following the renormalization group arguments of 
Schulz ||, an incommensurate umklapp term renormal- 
izes K Cl the charge field contribution to the exponents 
of correlation functions. The renormalized exponent, 
1/2 < K* < 1, tends to 1/2 for q — > q c and to 1 for 
<7> q'- 

We now address the question of deviation from per- 
fectly periodic impurity array. Since the spin part of 
the Hamiltonian is completely insensitive to the distance 
between impurities, the spin excitations will be gapped 
even in the case of a randomly spaced impurity array (i.e. 
disordered array). The charge field, on the other hand, 
is affected by the disorder. Following the method of Gia- 
marchi and Schulz || we conclude that a disordered ID 
Kondo array will be insulating. Therefore, the ID dis- 
ordered Kondo array will be a gapless charge insulator 
with a spin gap. 
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H c 
v F 



dx 



K + {d x <t> c - S) 2 



s 2 

dx—cos[f3(j) c ], 



where /? = y^r, 8 = f , u p = S 2 , £ = JTra 2 v F /J b x . 
Classically, this Hamiltonian is known to have an in- 
commensurate periodic soliton-lattice solution above 
a critical incommensurability q c 



2d J b ± 



TTVp ' 



(13) 



and to remain commensurate for q < q c . Thus, for q > q c 
there is a gapless Goldstone mode (acoustic vibrations of 
the soliton lattice). Note that the acoustic velocity u p 
is not equal to vp ; it gradually increases with 6. The 
solitons have a width L s — ^y^nvF / J±, so when the dis- 
tance between solitons become comparable with soliton 
width, there will be a "melting" crossover to the fully in- 
commensurate state discussed before. In our model this 
crossover will happen at incommensurability q' ~ 7t 2 <Zcj 
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